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Lignin  is  the  main  component  of  biomass  with  a  complex,  heterogeneous,  three-dimensional  polymeric 
structure  of  three  main  monolignols  (p-coumaryl,  coniferyl,  and  sinapyl  alcohol).  In  order  to  understand 
the  pyrolysis  mechanism  of  lignin  and  identify  the  chemical  pathways  for  the  formations  of  key  products 
during  pyrolysis,  the  pyrolysis  processes  of  /3-0-4  type  lignin  dimer  model  compound  1  (l-phenyl-2- 
phenoxy-1, 3-propanediol)  were  theoretically  investigated  by  employing  density  functional  theory  (DFT) 
methods  at  the  B3LYP/6-31G(d,p)  level.  Based  on  related  experimental  and  calculation  results  of  bond 
dissociation  energies  of  ft- 0-4  type  lignin  dimer,  three  possible  pyrolytic  pathways  (the  homolytic  cleav¬ 
age  of  C/j—O  bond,  the  homolytic  cleavage  of  Ca— bond  and  the  concerted  reactions)  were  proposed, 
the  activation  energies  of  each  reaction  step  were  calculated,  and  the  temperature  effect  on  pyrolysis  pro¬ 
cesses  was  analyzed.  The  calculation  results  indicate  that  the  homolytic  cleavage  reaction  of  C/j—O  bond 
and  concerted  reaction  pathways  (3)  could  be  the  major  reaction  channels,  and  the  homolytic  cleavage 
reaction  of  Ca—C/j  bond  and  concerted  reaction  pathways  (1)  and  (2)  could  be  the  competitive  reaction 
channels  in  pyrolysis  processes.  The  concerted  reactions  would  dominate  over  free-radical  homolytic 
reactions  at  lower  temperatures,  while  at  high  temperatures  the  free-radical  reaction  (C— 0  homolysis) 
would  dominate  over  the  concerted  reactions. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

With  increasing  concerns  about  energy  supply  and  environmen¬ 
tal  pollution  problems  caused  by  burning  fossil  fuels,  the  utilization 
of  renewable  biomass  resources  has  attracted  more  and  more 
attention.  The  sulfur,  nitrogen  and  ash  contents  in  biomass  are 
low  and  biomass  energy  is  a  renewable  and  clean  energy,  which 
can  reduce  the  greenhouse  effect,  making  the  biomass  energy 
study  a  hot  topic  worldwide  [1-8],  Biomass  liquefaction  technol¬ 
ogy  is  considered  to  be  the  most  promising  technology  in  biomass 
energy  utilizing  field.  Biomass  contains  varying  amounts  of  cellu¬ 
lose,  hemicellulose,  lignin  and  a  small  amount  of  extractives  [9], 
The  ratio  of  the  three  components  in  a  biomass,  which  is  closely 
related  to  biomass  conversion,  differs  greatly  along  various  sources 
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of  biomass  [10],  the  study  of  the  pyrolysis  mechanism  of  biomass 
components  is  therefore  essential. 

Lignin  is  the  second-most  abundant  naturally  occurring  biopoly¬ 
mer  and  a  by-product  of  the  pulping  process  [11],  Pyrolysis  of  lignin 
has  been  studied  by  some  people  over  the  past  several  decades, 
latridis  and  Gavalas  [12]  studied  the  pyrolysis  of  a  precipitated  kraft 
lignin  in  a  “captive  sample”  reactor  at  temperature  of  400-700  °C 
and  identified  a  few  compounds  by  gas  chromatography  including 
oxides  of  carbon,  hydrocarbons,  methanol,  acetone  and  single-ring 
phenols  (phenol,  guaiacol,  etc.).  Caballero  et  al.  [13]  used  a  semi¬ 
batch  flow  reactor  connected  in  series  between  a  Pyroprobe  1000 
and  a  gas  chromatographer  to  study  the  production  of  gases  from 
kraft  lignin  at  different  nominal  temperatures  (ranging  from  500 
to  900°C),  and  analyzed  yields  of  13  pyrolysis  products  (methane, 
ethane,  ethylene,  propylene,  propane,  hydrocarbons  with  four  car¬ 
bon  atoms,  benzene,  other  aromatics,  CO,  C02,  methanol  together 
with  formaldehyde,  acetaldehyde,  acetic  acid  and  water).  Jiang 
et  al.  [  1 4]  investigated  the  temperature  dependence  of  the  pyrolysis 
products  of  two  types  of  lignin  ( Alcell  lignin  and  Asian  lignin)  using 
PyGC-MS  and  identified  about  50  compounds  from  lignin  pyrolysis 
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at  a  temperature  range  of  400-800  °C.  The  phenolic  compounds 
yield  was  17.2  wt%  for  Alcell  lignin  and  individual  yield  of  most  of 
the  compounds  was  less  than  1  wt%.  Wang  et  al.  [15,16]  investi¬ 
gated  the  pyrolysis  characteristics  of  milled  wood  lignins  (MWLs) 
using  a  thermogravimetric  analyzer,  and  also  analyzed  the  FTIR 
profiles  of  main  pyrolysis  products  (C02,  CO,  methane,  methanol, 
phenols  and  formaldehyde). 

Lignin  model  compounds,  which  have  simple  structures  and 
product  distributions  compared  to  real  lignin,  have  been  studied 
in  order  to  understand  the  pyrolysis  mechanism  of  lignin.  Vuori 
and  Bredenberg  17]  studied  the  thermolyses  of  o-,  m-,  and  p- 
hydroxyanisoles  under  an  inert  atmosphere,  the  major  products 
formed  from  all  three  hydroxyanisoles  were  the  correspondingly 
substituted  dihydroxybenzenes  and  cresols,  and  the  reactivity  of 
the  lignin-related  o-hydroxyanisole  (guaiacol)  was  much  higher 
than  those  of  all  other  model  compounds  studied.  The  pyrolysis  of 
phenethyl  phenyl  ether  (PPE)  and  its  methoxy-substituted  deriva¬ 
tives  were  investigated  theoretically  by  Beste  et  al.  [18,19],  They 
found  that  the  initial  step  in  the  thermal  decomposition  of  PPE  is  the 
homolytic  cleavage  of  the  oxygen-carbon  bond,  and  carbon-carbon 
bond  cleavage  in  PPE  could  be  a  competitive  initial  reaction  under 
high-temperature  pyrolysis  conditions.  Chu  et  al.  [20]  investigated 
the  pyrolysis  behavior  of  a  fi-O-4  type  oligomeric  lignin  model  com¬ 
pound  at  a  temperature  range  from  250  to  550  °C  and  proposed  a 
free  radical  reaction  pathway  to  explain  the  product  chemistry.  To 
obtain  information  on  the  thermal  reactivities  of  lignin  aromatic 
nuclei,  guaiacyl  and  syringyl  types,  thermal  reactions  of  guaia¬ 
col  (2-methoxyphenol)  and  syringol  (2,6-dimethoxyphenol)  were 
compared  in  a  closed  ampoule  reactor  (N2/400-600°C/40-600s) 
by  Asmadi  et  al.  [21],  Their  experimental  results  indicated  that 
the  O-CH3  bond  homolysis  occurring  at  temperatures  higher  than 
400  C  was  the  rate-determining  step,  and  the  radical  induced 
rearrangement  (ipso-substitution),  demethoxylation  and  radical 
coupling  reactions  to  form  higher  MW  products  and  methylated 
aromatics  were  the  subsequent  reactions. 

Lignin  is  a  complex,  heterogeneous,  three-dimensional  polymer 
and  naturally  produced  by  the  enzymatic  polymerization  of  three 
main  monolignols  (p-coumaryl,  coniferyl,  and  sinapyl  alcohol), 
which  are  interconnected  by  a  variety  of  linkages  (j3- 0-4,  a- 0-4, 
4-0-5,  5-5,  jS-1,  etc.).  The  fi- 0-4  bond  is  the  major  type  of  linkage 
which  occupies  48-60%  of  the  total  linkages  depending  on  the  type 
of  wood  [22],  Unlike  the  pyrolysis  of  cellulose  and  hemicelluloses 
[23-25],  the  complexity  of  the  lignin  structure  and  its  high  molec¬ 
ular  weight  present  make  the  lignin  pyrolysis  more  complex,  and 
it  is  difficult  to  analyze  detailed  mechanism  through  experiments. 
So  people  gradually  pay  more  attention  to  employing  theoretical 
methods  to  study  the  pyrolysis  mechanism  and  to  forecast  possible 
reaction  pathways  [26-29],  In  order  to  further  improve  the  under¬ 
standing  of  detailed  mechanism  of  lignin  pyrolysis  and  formation 
mechanism  of  its  main  products,  the  pyrolysis  processes  of  fi-O-4 
type  lignin  dimer  model  compound  1  (l-phenyl-2-phenoxy-l,3- 
propanediol,  shown  in  Fig.  1 )  have  been  theoretically  investigated 
by  using  density  functional  theory  methods  at  B3LYP/6-31G(d,p) 
level  in  this  paper. 


2.  Calculation  methods 

The  geometries  of  reactant,  intermediates,  transition  states 
and  products  have  been  optimized  by  employing  the  hybrid  den¬ 
sity  functional  B3LYP  (Becke’s  three  parameter  gradient  corrected 
exchange  functional  [30]  with  the  gradient  corrected  correlation 
functional  of  Lee  et  al.  [31])  with  the  6-31G(d,p)  basis  set.  Station¬ 
ary  points  on  the  potential  energy  surface  of  the  reacting  system 
were  fully  optimized  using  gradient  minimization  techniques,  fol¬ 
lowed  by  evaluating  harmonic  vibration  frequencies  to  characterize 


their  nature  as  minima  (all  real  frequencies)  or  first-order  saddle 
points  (namely  transition  states,  only  one  imaginary  frequency).  No 
constrain  in  the  initial  molecular  geometry  for  searching  symmet¬ 
ric  and  asymmetric  transition  states  was  used  in  their  optimization 
procedure.  Intrinsic  reaction  coordinate  (IRC)  [32]  calculations  have 
been  performed  to  identify  the  minimum  energy  path  from  a  tran¬ 
sition  state  to  the  two  corresponding  local  minima.  Activation 
energies  (the  reaction  energy  barriers)  for  reactions  were  estimated 
as  the  relative  energies,  including  zero-point  energy  correction 
(ZPE),  between  the  transition  state  and  the  reactant.  The  standard 
thermodynamic  changes  were  energy  differences  between  the 
reactants  and  the  products,  including  the  ZPE.  For  free-radical  reac¬ 
tions,  the  bond  dissociation  energies  may  serve  as  approximations 
to  the  activation  energies.  Gaussian  03  suites  of  programs  were 
employed  for  all  calculations  [33], 

3.  Results  and  discussion 

3.2.  Bond  dissociation  energies  for  fi-O-4  type  lignin  dimer  model 
compound 

Bond  dissociation  energies  ( A£)  are  fundamental  to  understand 
many  chemical  processes  because  they  exemplify  the  reaction 
energy  of  a  homolytic  cleavage  of  the  considered  bond.  The 
homolytic  bond  dissociation  is  an  important  reaction  step  in  lignin 
pyrolysis.  Fig.  1  shows  the  homolytic  cleavage  ways  of  fl-O-4  type 
lignin  dimer  model  compound  and  their  bond  dissociation  ener¬ 
gies  calculated  at  B3LYP/6-31G(d,p)  level.  The  bond  dissociation 
energy  of  Cp—O,  245.3  kj/mol,  is  the  lowest,  and  that  of  Ca—Cp, 
259.2  kj/mol,  is  the  second  lowest.  The  order  of  A E  is  as  follows: 
C/3—O  <  Ca-Cp  <  C„-OH  <  Cp-Cy  <  Caromatic-C„<  Caromatic-0,  there¬ 
fore,  the  initial  step  in  the  thermal  decomposition  of  fi- 0-4  type 
lignin  dimer  model  compound  could  be  the  homolytic  cleavage  of 
the  Cp—O  bond,  and  Ca—Cp  bond  cleavage  could  be  a  competitive 
initial  reaction. 

3.2.  Pyrolysis  mechanism  of  fi-O-4  type  lignin  dimer  model 
compound 

Lignin  pyrolysis  experimental  results  [14,34,35]  indicated  that 
the  main  pyrolysis  products  were  single-ring  phenolic  compounds, 
methanol,  formaldehyde,  acetaldehyde,  CO,  CH4,  etc.  Britt  et  al. 
[22]  studied  the  fast  vacuum  pyrolysis  of  methoxy-substituted  fi- 
O-4  lignin  model  compounds  and  proposed  a  complex  reaction 
pathway  which  was  dominated  by  free-radical  reactions,  molec¬ 
ular  rearrangements,  and  concerted  elimination  reactions.  In  this 
paper,  we  propose  three  possible  pyrolysis  ways:  the  homolytic 
cleavage  of  Cp—O  bond,  the  homolytic  cleavage  of  Ca—Cp  bond  and 
the  concerted  reactions. 

3.2.1.  The  homolytic  cleavage  ofCp—O  bond 

According  to  the  above  calculation  results  of  bond  dissociation 
enthalpies  of  fi-O-4  type  lignin  dimer  model  compound,  the  bond 
dissociation  energy  of  Cp—O  bond  is  the  lowest,  245.3  kj/mol,  and 
the  homolytic  cleavage  of  Cp—O  bond  is  the  major  pyrolysis  way 
in  pyrolysis  of  lignin.  Fig.  2  gives  the  possible  reaction  pathways 
for  the  homolytic  cleavage  of  Cp—O  bond,  the  potential  energy  pro¬ 
file  along  reaction  pathways  is  shown  in  Fig.  3,  and  the  optimized 
geometries  of  reactant,  intermediates,  transition  states,  and  prod¬ 
ucts  in  pyrolysis  processes  of  the  homolytic  cleavage  of  C^— 0  bond 
are  shown  in  Fig.  4. 

Model  compound  1  can  decompose  into  radicals  2  and  3,  then 
the  radical  3  is  converted  into  phenol  4  through  hydrogenation 
reaction,  which  is  exothermic  with  an  energy  of-340.3  kj/mol.  Rad¬ 
ical  2  can  undergo  a  dehydrogenation  to  form  phenylallyl  alcohol 
5  via  the  transition  state  T1  with  an  energy  barrier  of  143.3  kj/mol, 
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Fig.  1.  Homolytic  cleavage  ways  of  fi-O-4  type  lignin  dimer  model  compound  and  their  bond  dissociation  enthalpies  ( A£). 


absorbing  an  energy  of  124.2  kj/mol,  or  it  can  undergo  a  hydro¬ 
genation  to  form  phenyl  propanediol  9,  releasing  an  energy  of 
378.4  kj/mol.  Compound  5  can  tautomerize  to  form  phenyl  acetone 
alcohol  6  via  T2  with  an  activation  energy  of  200.1  kj/mol  through 


intramolecular  hydrogen  transfer.  The  tautomerization  reaction  is 
exothermic  and  compound  6  is  more  stable  than  its  enol  form 
compound  5.  Compound  6  can  decompose  into  acetophenone  7 
and  formaldehyde  via  transition  state  T3  with  an  energy  barrier 


Fig.  2.  The  proposed  reaction  pathways  in  pyrolysis  processes  of  the  homolytic  cleavage  of  Cp— 0  bond. 
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Fig.  3.  Potential  energy  profile  along  reaction  pathways  in  pyrolysis  processes  of  the  homolytic  cleavage  of  Cp— O  bond. 


of  345.5  kj/mol,  or  it  can  undergo  a  dehydration  reaction  to  form 
phenyl  acrylketone  8  and  H2O  via  transition  state  T4  with  an  energy 
barrier  of  241.7  kj/mol.  Phenyl  propanediol  9  can  decompose  via 
three  pathways.  In  pathway  (1 ),  compound  9  decomposes  to  form 
phenethyl  alcohol  10  and  formaldehyde  via  transition  state  T5  with 
an  energy  barrier  of  359.3  kj/mol.  In  pathway  (2),  compound  9 
decomposes  to  form  benzaldehyde  11  and  ethyl  alcohol  12  via  tran¬ 
sition  state  T6  with  an  energy  barrier  of  324.2  kj/mol.  In  pathway 
(3),  compound  9  undergoes  a  dehydration  reaction  to  form  com¬ 
pound  13  and  H2O  via  transition  state  T7  with  an  energy  barrier  of 
260.2  kj/mol,  and  compound  13  tautomerizes  to  form  compound 
14  via  T8  with  an  energy  barrier  of  272.5  kj/mol  through  intra¬ 
molecular  hydrogen  transfer.  From  the  optimized  structure  of  T8 
in  Fig.  4,  there  is  a  transfer  of  hydrogen  atom  from  CY  to  Q,,  mean¬ 
while  Ca—Cfj  double  bond  lengthens  to  transform  into  single  bond 
and  Cjj—CY  bond  shortens  to  transform  into  double  bond.  Com¬ 
pound  14  can  further  tautomerize  to  form  benzenepropanal  15  via 
T9  with  an  activation  energy  of  231.0  kj/mol.  Benzenepropanal  15 
can  decarbonylate  to  form  phenylethane  16  and  CO  via  T10  with 
an  energy  barrier  of  342.2  kj/mol. 

Through  comparisons  of  energy  barriers  of  each  reaction  step 
in  pyrolysis  processes  of  the  homolytic  cleavage  of  Cjg— O  bond,  we 
would  infer  that  the  major  products  of  model  compound  1  pyrolysis 
appear  to  be  phenol  4,  phenyl  acetone  alcohol  6,  phenyl  acrylketone 
8,  phenyl  propanediol  9,  benzenepropanal  15,  etc. 

3.2.2.  The  homolytic  cleavage  ofCa—Cp  bond 

Fig.  5  gives  the  possible  reaction  pathways  for  the  homolytic 
cleavage  of  Ca—Cp  bond,  the  potential  energy  profile  along  reac¬ 
tion  pathways  is  shown  in  Fig.  6,  and  the  optimized  geometries  of 
reactants,  intermediates,  transition  states,  and  products  in  pyroly¬ 
sis  processes  of  the  homolytic  cleavage  of  Ca—Cp  bond  are  shown  in 
Fig.  7.  The  radicals  17  and  18  can  be  formed  through  the  homolytic 
cleavage  of  the  Ca—Cp  bond.  The  radical  18  can  be  converted  into 
phenylcarbinol  19  through  hydrogenation  reaction,  releasing  an 
energy  of  332. 5  kj/mol,  or  into  benzaldehyde  11  through  dehydro¬ 
genation  reaction  via  transition  state  Til  with  an  energy  barrier 
of  157.61<J/mol.  The  radical  17  can  decompose  through  three  path¬ 
ways.  In  pathway  (1),  radical  17  loses  a  hydrogen  atom  through 
/3-scission  to  form  phenoxy  vinyl  alcohol  25  via  transition  state 
T12  with  an  energy  barrier  of  217.6kJ/mol.  Phenoxy  vinyl  alcohol 
25  tautomerizes  to  form  phenoxy  acetaldehyde  28  via  T13  with 


an  energy  barrier  of  225.9  kj/mol  through  intramolecular  hydro¬ 
gen  shift.  Phenoxy  acetaldehyde  28  may  decarbonylate  to  produce 
anisole  21  via  T14  with  a  high  energy  barrier  of  340.6  kj/mol. 
In  pathway  (2),  radical  17  decomposes  to  form  radical  24  and 
formaldehyde  via  T15  with  a  low  energy  barrier  of  83.71<J/mol. 
Radical  24  tautomerizes  to  form  radical  26  via  T16  with  an  energy 
barrier  of  86.2  kj/mol,  then  the  radical  26  is  further  converted  into 
radical  27  through  the  cleavage  of  the  C— O  bond  and  releases  an 
energy  of  51.1  kj/mol.  The  radical  27  can  pick  up  a  hydrogen  atom 
to  become  phenylcarbinol  19,  releasing  an  energy  of  404.2  kj/mol, 
or  it  can  lose  a  hydrogen  atom  to  form  benzaldehyde  11  via  transi¬ 
tion  state  T17  with  a  low  energy  barrier  of  80.9  kj/mol.  In  pathway 
(3),  radical  17  undergoes  a  hydrogenation  to  form  phenoxy  alcohol 
20,  releasing  an  energy  of  383.1  kj/mol.  Phenoxy  alcohol  20  may 
decompose  to  form  anisole  21  and  formaldehyde  via  T18  with  a 
high  energy  barrier  of  367.7  kj/mol,  or  it  may  form  phenol  4  and 
vinyl  alcohol  22  via  T19  with  an  energy  barrier  of  257.1  kj/mol, 
which  is  similar  to  a  dehydration  reaction.  From  the  optimized 
structure  of  T19  in  Fig.  7,  C^— O  and  Cy— H  bonds  rupture,  mean¬ 
while  Cp—CY  bond  shortens  to  transform  into  double  bond.  Vinyl 
alcohol  22  can  tautomerize  to  form  acetaldehyde  23  via  T20  with 
an  energy  barrier  of  224.0  kj/mol  through  intramolecular  hydrogen 
transfer  and  release  an  energy  of  50.8  kj/mol. 

Through  comparisons  of  energy  barriers  of  each  reaction  step  in 
pyrolysis  processes  of  the  homolytic  cleavage  of  Ca— bond,  we 
would  infer  that  the  major  products  of  model  compound  1  pyrol¬ 
ysis  appear  to  be  phenol  4,  benzaldehyde  11,  phenylcarbinol  19, 
phenoxy  alcohol  20,  acetaldehyde  23,  phenoxy  acetaldehyde  28, 
etc. 

3.2.3.  The  concerted  reactions 

Fig.  8  gives  four  possible  concerted  reaction  pathways  in  pyrol¬ 
ysis  processes  of  /1-0-4  type  lignin  dimer  model  compound,  the 
potential  energy  profile  along  reaction  pathways  is  shown  in  Fig.  9, 
and  the  optimized  geometries  of  reactant,  intermediates,  transi¬ 
tion  states,  and  products  in  concerted  reaction  pathways  are  shown 
in  Fig.  10.  In  concerted  reaction  pathway  (1),  model  compound 
1  decomposes  directly  to  form  phenol  4  and  intermediate  29  via 
transition  state  T21  with  an  energy  barrier  of  257.1  kj/mol  through 
concerted  elimination  reaction  without  the  intermediacy  of  free 
radicals.  From  the  optimized  structure  ofT21  in  Fig.  10,  bond  Cp—O 
and  bond  O— H  on  C„  rupture,  meanwhile  phenol  4  is  formed. 
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Fig.  4.  Optimized  geometries  of  reactant,  intermediates,  transition  states,  and  products  in  pyrolysis  processes  of  the  homolytic  cleavage  of  C^— 0  bond  (unit:  A). 
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Fig.  6.  Potential  energy  profile  along  reaction  pathways  in  pyrolysis  processes  of  the  homolytic  cleavage  of  Ca—Cp  bond. 
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Intermediate  29  tautomerizes  to  form  phenyl  acetone  alcohol  6  via 
T25  with  an  energy  barrier  of  258.1  kj/mol  through  intramolecular 
hydrogen  transfer,  releasing  an  energy  of  111.6kJ/mol.  From  the 
optimized  structure  of  T25  in  Fig.  10,  there  is  a  transfer  of  hydro¬ 
gen  atom  from  Ca  to  Cp,  meanwhile  O—Cp  bond  ruptures  and 
0— C„  bond  shortens  to  transform  into  double  bond.  Compound  6 
may  further  decompose  into  acetophenone  7  and  formaldehyde 
via  transition  state  T3  with  a  high  energy  barrier  of  345.5  kj/mol. 
In  concerted  reaction  pathway  (2),  phenol  4  and  3-hydroxy-3- 
phenyl-propaldehyde  30  are  directly  formed  via  transition  state 
T22  with  an  energy  barrier  of  240.1  kj/mol  through  concerted 
elimination  reaction  without  the  intermediacy  of  free  radicals.  It 
can  be  seen  from  the  optimized  structure  of  T22  in  Fig.  10  that 
Cp—O  bond  and  O— H  bond  on  Cy  rupture,  meanwhile  0—Cy  bond 
shortens  to  transform  into  double  bond  and  phenol  4  is  formed. 
3-hydroxy-3-phenyl-propaldehyde  30  may  decarbonylate  to  form 


phenethyl  alcohol  10  and  CO  via  T26  with  a  high  activation  energy 
of  335.0  kj/mol,  and  phenethyl  alcohol  10  may  undergo  a  demeth- 
ylation  to  form  benzaldehyde  11  and  CH4  via  transition  state  T27 
with  an  energy  barrier  of  325.8  kj/mol  through  concerted  elimina¬ 
tion  reaction.  In  concerted  reaction  pathway  (3),  model  compound 
1  decomposes  directly  to  form  phenol  4  and  phenylallyl  alcohol  5 
via  transition  state  T23  with  an  energy  barrier  of  21 9.3  kj/mol.  It  can 
be  seen  from  the  optimized  structure  ofT23  in  Fig.  10  that  Cp—O 
bond  and  O— H  bond  on  Cff  rupture,  meanwhile  Ca—Cp  bond  has 
been  shortened  from  1.546  A  to  1.421  A  to  transform  into  double 
bond  and  phenol  4  and  phenylallyl  alcohol  5  are  formed.  Phenylallyl 
alcohol  5  can  tautomerize  to  form  phenyl  acetone  alcohol  6  via  T2 
with  an  activation  energy  of  200.1  kj/mol  through  intramolecular 
hydrogen  transfer,  releasing  an  energy  of  55.91<J/mol.  Phenyl  ace¬ 
tone  alcohol  6  may  further  decompose  into  acetophenone  7  and 
formaldehyde.  In  concerted  reaction  pathway  (4),  benzaldehyde 
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Fig.  8.  Four  possible  concerted  reaction  pathways  in  pyrolysis  processes  of  /J-O-4  type  lignin  dimer  model  compound. 
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Fig.  9.  Potential  energy  profile  along  reaction  pathways  in  four  concerted  reaction  pathways. 
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11  and  phenoxy  alcohol  20  are  directly  formed  via  transition  state 
T24  with  a  high  energy  barrier  of  343.9  kj/mol  through  concerted 
elimination  reaction.  It  can  be  seen  from  the  optimized  structure  of 
T24  in  Fig.  10  that  Ca—Cp  bond  and  O— H  bond  on  Ca  rupture,  mean¬ 
while  0— Cff  bond  shortens  to  transform  into  double  bond.  Phenoxy 
alcohol  20  may  further  decompose  through  the  proposed  reaction 
pathways  in  Fig.  5. 

Comparing  energy  barriers  of  initial  reaction  steps  in  four  pos¬ 
sible  concerted  reaction  pathways,  it  is  easy  to  know  that  energy 
barriers  of  initial  reaction  steps  of  concerted  reaction  pathways 

(1) ,  (2)  and  (3)  are  obviously  lower  than  that  of  concerted  reac¬ 
tion  pathway  (4),  so  concerted  reaction  pathways  (1),  (2)  and  (3) 
are  kinetically  favorable,  and  concerted  reaction  pathway  (4)  is 
unlikely  to  occur  in  pyrolysis  processes  of ) 3-0-4  type  lignin  dimer 
model  compound.  For  concerted  reaction  pathways  (1),  (2)  and 
(3),  the  energy  barrier  of  initial  reaction  step  of  concerted  reaction 
pathway  (3)  is  lower  than  those  of  concerted  reaction  pathways  (1 ) 
and  (2),  consequently  concerted  reaction  pathway  (3)  could  be  the 
major  reaction  channel  and  concerted  reaction  pathways  ( 1 )  and  (2 ) 
could  be  the  competitive  reaction  channels  in  pyrolysis  processes 
of  model  compound  1. 

Through  comparisons  of  energy  barriers  of  initial  reaction  steps 
in  three  possible  pyrolysis  ways:  the  homolytic  cleavage  reaction 
of  Cp—O  bond,  the  homolytic  cleavage  reaction  of  Ca—Cp  bond  and 
the  concerted  reactions,  it  can  be  seen  that  the  homolytic  cleav¬ 
age  reaction  of  Cp—O  bond  and  concerted  reaction  pathways  (3) 
could  be  the  major  reaction  channels,  and  the  homolytic  cleavage 
reaction  of  Ca—Cp  bond  and  concerted  reaction  pathways  (1)  and 

(2)  could  be  the  competitive  reaction  channels  in  pyrolysis  pro¬ 
cesses  of  /3-0-4  type  lignin  dimer  model  compound  1.  Hence,  the 
main  products  of  model  compound  1  pyrolysis  could  be  pheno¬ 
lic  compounds  such  as  phenol  4,  phenyl  acetone  alcohol  6,  phenyl 
acrylketone  8,  phenyl  propanediol  9  and  benzenepropanal  15,  and 


the  main  competitive  products  could  be  benzaldehyde  11,  phenyl- 
carbinol  19,  phenyl  ether  such  as  phenoxy  alcohol  20,  phenoxy 
vinyl  alcohol  25  and  phenoxy  acetaldehyde  28.  The  above  analy¬ 
sis  results  are  in  accordance  with  the  related  experimental  results 

[14,34,35], 

3.3.  Temperature  effect  on  pyrolysis  processes  of  / 3-0-4  type 
lignin  dimer  model  compound 

In  order  to  understand  the  temperature  effect  on  pyrolysis  pro¬ 
cesses  of  (6-0-4  type  lignin  dimer  model  compound  1,  activation 
energies  of  initial  reaction  steps  in  pyrolysis  processes  of  model 
compound  1  at  different  temperatures  (298,  450,  600,  750,  900, 
1100,  1300,  1600,  1900  K)  were  calculated  at  B3LYP/6-31G(d,p) 
level.  Table  1  summarizes  the  activation  energies  of  initial  reac¬ 
tion  steps  at  different  temperatures,  and  the  relationship  between 


Table  1 

Activation  energies  of  initial  reaction  steps  in  pyrolysis  processes  of  model  com¬ 
pound  1  at  different  temperatures. 


Temperature 

(K) 

Activation  energies  (kj/mol) 

Reaction  Reaction 

step  1  step  2 

Reaction 
step  3 

Reaction 
step  4 

Reaction 
step  5 

298 

245.3 

259.2 

257.1 

240.1 

219.3 

450 

244.2 

258.9 

256.9 

239.9 

220.1 

600 

242.4 

257.6 

256.6 

239.6 

220.5 

750 

240.1 

255.6 

256.2 

239.3 

220.7 

900 

237.5 

253.1 

255.7 

239.1 

220.8 

1100 

233.7 

249.3 

254.8 

238.9 

220.6 

1300 

229.5 

245.1 

253.8 

238.5 

220.1 

1600 

223.0 

238.4 

252.3 

237.8 

219.0 

1900 

216.2 

231.4 

250.5 

236.8 

218.0 

Reaction  step  1 : 

1^2  +  3; 

reaction  step  2: 

1  17  +  18; 

reaction  step  3:  1  ^T21; 

reaction  step  4: 1  ->  T22;  reaction  step  5:1  >  T23. 
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Fig.  11.  Relationship  between  activation  energies  of  initial  reaction  steps  in  pyrol¬ 
ysis  processes  of  model  compound  1  and  temperature  (reaction  step  1:  1^2  +  3; 
reaction  step  2:  1  >17+18;  reaction  step  3:  1  >T21;  reaction  step  4:  1  ->  T22; 
reaction  step  5: 1  -> T23). 


(2)  The  homolytic  cleavage  reaction  of  Cp—O  bond  and  concerted 
reaction  pathways  (3 )  could  be  the  major  reaction  channels,  and 
the  homolytic  cleavage  reaction  of  Ca—Cp  bond  and  concerted 
reaction  pathways  ( 1 )  and  (2)  could  be  the  competitive  reaction 
channels  in  pyrolysis  processes. 

(3)  The  main  products  of  model  compound  1  pyrolysis  could  be 
phenolic  compounds  such  as  phenol  4,  phenyl  acetone  alcohol 
6,  phenyl  acrylketone  8,  phenyl  propanediol  9  and  benzene- 
propanal  15,  and  the  main  competitive  products  could  be 
benzaldehyde  11,  phenylcarbinol  19,  phenyl  ether  such  as 
phenoxy  alcohol  20,  phenoxy  vinyl  alcohol  25  and  phenoxy 
acetaldehyde  28. 

(4)  The  concerted  reactions  would  dominate  over  free-radical 
homolytic  reactions  at  lower  temperatures,  while  at  high  tem¬ 
peratures  the  free-radical  reaction  (C— 0  homolysis)  would 
dominate  over  the  concerted  reactions. 

Acknowledgments 


activation  energy  and  temperature  is  depicted  in  Fig.  11.  It  is 
seen  from  Fig.  1 1  that  activation  energies  of  free-radical  reactions 
(steps  1  and  2)  are  obviously  temperature-dependent,  but  acti¬ 
vation  energies  of  concerted  reactions  (steps  3-5)  show  limited 
temperature-dependency.  With  the  increase  of  temperature,  the 
activation  energies  of  free-radical  reactions  ( steps  1  and  2 )  decrease 
obviously,  while  the  activation  energies  of  concerted  reactions 
(steps  3-5)  show  little  change.  Jarvis  et  al.  [36;  studied  the  pyroly¬ 
sis  of  PPE  as  lignin  model  compound  in  a  wide  temperature  range 
of  300-1350  °C  to  determine  the  importance  of  concerted  reaction 
and  free-radical  reaction  pathways,  and  experimental  results  indi¬ 
cate  that  the  free-radical  reaction  (C— O  homolysis)  is  significant  at 
high  temperatures  (>1000  °C),  whereas  the  concerted  reactions  are 
significant  at  lower  temperatures,  and  the  C— C  homolysis  mech¬ 
anism  is  minor  at  all  temperatures.  The  structure  of  /3-0-4  type 
lignin  dimer  model  compound  1  is  similar  to  PPE,  but  more  com¬ 
plex  than  PPE.  We  can  see  from  Table  1  that  at  293  K  the  activation 
energy  of  concerted  reaction  (step  5)  is  the  lowest,  however,  with 
the  increase  of  temperature,  the  activation  energy  of  free-radical 
reaction  (C— O  homolysis )  becomes  the  lowest  at  1 900  K.  Our  calcu¬ 
lation  results  indicate  that  the  concerted  reactions  would  dominate 
over  free-radical  homolytic  reactions  at  lower  temperatures,  while 
at  high  temperatures  the  free-radical  reaction  (C— O  homolysis) 
would  dominate  over  the  concerted  reactions,  which  is  basically 
consistent  with  the  result  of  the  literature  [36], 

4.  Conclusions 

In  this  work,  the  pyrolysis  mechanisms  of  /3-0-4  type  lignin 
dimer  model  compound  1  have  been  theoretically  investigated  by 
employing  DFT  methods  at  B3LYP/6-31G(d,p)  level.  Three  possi¬ 
ble  pyrolytic  pathways  (the  homolytic  cleavage  of  Cg— O  bond,  the 
homolytic  cleavage  of  Ca—Cp  bond  and  the  concerted  reactions) 
were  proposed  and  the  geometries  of  reactants,  transition  states, 
intermediates  and  products  were  fully  optimized.  The  activation 
energies  of  each  reaction  step  were  calculated  and  the  temperature 
effect  on  pyrolysis  processes  was  analyzed.  The  main  conclusions 
are  made  as  follows: 

(1)  The  order  of  magnitude  of  bond  dissociation  energies  for 
model  compound  1  is  as  follows:  C^— 0<Cq.— C^<Cq.— OH< 
C^-Cy  <  Caromatic— Cff<  Caromatic— o,  so  the  initial  step  in  the 
thermal  decomposition  of  model  compound  1  could  be  the 
homolytic  cleavage  of  the  Cp— O  bond,  and  Ca—Cp  bond  cleavage 
could  be  a  competitive  initial  reaction. 
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